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There is concern that genetically modified (GM) plants may have adverse affects on the arthropod 
biodiversity comprising agricultural landscapes. The present study report on a two year field experimental 
test of whether four different genotypic lines, some are novel with no previous field tests, of GM maize 
hybrids alter the structure of arthropod food webs that they harbour, relative to non-GM maize (control) 
that is widely used in agriculture. The different GM genotypes produced either Bt toxins, conferred 
glyphosate tolerance or a combination of the two traits. Quantitative food web analysis, based on short-term 
assessment assigning a total of 243,896 arthropod individuals collected from the treatments to their 
positions in food webs, revealed that complex and stable food webs persisted in each maize treatment. 
Moreover, food web structure remained relatively unchanged by the GM-genotype. The results suggest that 
at least in short-term period these particular GM maize genotypes will not have adverse effects on arthropod 
biota of agricultural landscapes. 



There is growing appreciation that the genotype of plants can be an important determinant of associated 
arthropod food webs' These genetic effects on arthropod food webs may also have cascading effects on the 
functioning of ecosystems''''' ''. Given such relationships, it stands to reason that the release of genetically 
modified (GM) crops into the environment could disrupt arthropod food web structure and ecosystem function- 
ing, especially through impacts on non-target arthropod species''-'. However, information on any potential link 
between GM genotypes and arthropod food web structure is scarce if nonexistent, thereby precluding the 
possibility to better inform policy about the broad impacts of GM crops are on the environment. The present 
study reports on the effects of four different genotypic lines of GM maize from which two with extra glyphosate 
treatment hybrids relative to non-GM maize on the structure of arthropod food webs that they harbour (Table 1). 
Through quantitative food web analysis the study offer insight into how food web properties such as bottom-up 
(Kbu.i) and top-down (K,;j.i) indexes and connectance are affected by the GM maize genotype. 

The approach is to examine the potential for tissue or genetic material from GM maize to enter no-target 
species via food chain effects or whether GM lines would disrupt food chain structure to the point where 
transmission is altered, in light of ecological observations that plant genotypes differentially affect arthropod 
species community structure. A potentially important route of exposure is through transgenic proteins that are 
bioaccumulated up trophic levels in a food chain after they are consumed by herbivores feeding on maize" '" that 
are then consumed by their predators (or consumption of pollen if prey are scarce". Such bioaccumulation could 
have deleterious effects on predator survival, longevity and development, thereby shortening food chains on GM 
plants relative to non-GM plants, depending on genotype-specific toxicity. The presence of such affects have been 
examined in several laboratory studies with a host of arthropod species and Bt protein from crops'" '"' '■*. These 
generally shown that effects of Bt protein may be negligible. But, these laboratory studies focus on single predator 
and prey species, and thus may fail to measure any potential sub-lethal effects that only become manifest when 
species are juxtaposed and interact with other species within their natural food webs, as is the case when 
examining effects of pesticides on predator-prey performance between lab and field settings'"""". Field studies 
are therefore needed to comprehensively assess the effects of Bt protein exposure to many insect species. 
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More than 50 field studies have been conducted in commercial and 
experimental fields to evaluate the impact of GM crops that contain 
transgen from Bacillus thuringiensis (Bt crops) on beneficial non- 
target organisms'^. The vast majority of studies have focused on Bt 
maize commercially grown since 1996" '^. These studies revealed 
minor, transient and inconsistent effects on non-target organisms 
when compared with non-Bt controls' '^. All important phytophag- 
ous and predators species were previously considered non-target 
organisms on which possible harmful effects of Bt toxins may have 
been directly (phytophagous) or indirectly (predators) detected" 
Currently, there are no evidence to prove that Bt proteins accumulate 
in prey tissues'". Other studies individually analyzing these groups 
also revealed no deleterious effects of Bt maize on mortality (studies 
on lacewing larvae and lady beetles'^, adult and juvenile Theridion 
impressum spiders'"'^, field and laboratory studies on Orius 
spp., lacewings and Stethorus punctUlum', staphylinidae^" and 
Lepidoptera^'. There were also no significant changes in heteropteran 
predator densities exposed with herbicide tolerant maize varieties^^. 
Field studies can indeed assess the effects of exposure on many spe- 
cies^''^"*, but quantitative food web analyses of arthropods on GM 
maize crops can provide a powerful tool to explore the structure of 
ecological communities and their responses to environmental factors 
that may not be revealed by conventional analyses of species presence 
or absence^^. During this two-year field study detailed arthropod 
assessments were made on GM maize treatments and on non-GM 
controls (Table 1) to test the hypothesis that a similar and stable food- 
web exists on all maize crops. Although the total number of arthro- 
pods used for analyses exceeded 240,000 these can only be considered 
as a short-term (or acute) impact of transgenic crops on arthropod 
food webs. Therefore to increase the sensibility of the methods food- 
web analyses were mainly focus on a subset of organisms that are: 

1. Most abundant during the maize vegetation period (but before 
poUen spreading to exclude errors by cross-poUination). 

2. Their cumulative abundance as trophic groups (herbivores, pre- 
dators) for each GM and non-GM crops exceeded 95%. 

3. Were identified as test organisms in GM crops by EFSA^''. 

The two target organisms western corn rootworm {Diabrotica 
virgifera virgiferd) and the European corn borer {Ostrinia nubilalis) 
were also included in food-web analyses. 

Study system. GM maize is grown annually on more than 35 million 
hectares globally^'. Among EU countries, maize is cultivated on 13 
million hectares, which represents 13% of the cultivated area in the 
EU or 8% of maize production area worldwide^". Currently, only one 



GM maize, MON810 that expresses active Bt toxin that protects it 
against insect pests, is legally grown in open field cultivation for 
commercial purposes in Europe. 

The most important insect pests of non-GM field maize in Europe 
are the western corn rootworm {Diabrotica virgifera virgifera) and 
the European corn borer {Ostrinia nubilalis). Other insects, such as 
aphids and thrips, occur in all maize crops and may occasionally 
cause losses. The most important predators of these pests and other 
arthropod herbivores are lady beetles, predatory thrips, rove beetles 
(Staphylinidae), spiders, Syrphid larvae and minute pirate bugs 
(Orius spp)^". 

According to the ISSAA^', the main applications for GM maize at 
present are: (1) pest resistance; (2) herbicide tolerance; and (3) a 
combination of the two. Pest resistant maize contains toxins against 
the western corn rootworm and the European corn borer. Herbicide 
tolerant maize, especially glyphosate-resistance genotypes are 
important for early growth stages of maize, because application of 
this broad spectrum herbicide reduces competition maize plants and 
weeds, and had no negative effects on maize plants. 

The GM varieties studied in the present experiments contain tox- 
ins against insect pests, or proteins that confer tolerance to glypho- 
sate, or in combination confer tolerance against both insects and 
glyphosate. Four GM varieties whose protein make-up conferred 
different combinations of resistance to pest species and tolerance 
to glyphosate were evaluated (Table 1). Resistance in different GM 
varieties was conferred through different Bacillus thuringiensis 
insecticidal crystal (Cry) proteins in different forms that can specif- 
ically target different insect pest species or the combinations of those 
pests. Some GM varieties also contain proteins that confer tolerance 
to glyphosate herbicides in addition to insecticidal proteins (Table 1). 

Results 

The total number of individuals used for food-web construction 
varied between treatments, ranging from the highest with 32,237 
in Lepidoptera resistant and glyphosate tolerant + glyphosate 
treated blocks and the lowest with 24,972 in Coleoptera and 
Lepidoptera resistant blocks (Table 2). Different treatments had very 
similar food web structure with relative constant numbers of trophic 
species (S) and trophic links (L) between species (Figures 1 A-H). 
Predatory groups of all treatments were dominated by ladybird bee- 
tles, Lacewing larvae, predatory thrips, rove beetles {Staphylinidae), 
minute pirate bugs {Orius), dipteral syrphidae, and spiders. 
Phytophagous groups were dominated by aphids, Cicadellidae, 
thrips, flea beetles, and Diabrotica virgifera virgifera was also present. 
The most frequent weed species were Echinochloa crus-galli. 
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Table 2 | The number of total 


individuals used for the quantitative food-web constructions 
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Cynodon dactylon. Sorghum halapense. Convolvulus arvensis, 
Amaranthus retroflexus, and Rubus caesius. (Figures 1 A-H). All 
these non-targett organisms were defined as test organisms in GM 
crops. Specialist insects that depend more on target pests may be the 
exception that Bt maize does not impact non-target organisms. This 
is particularly true for parasitoids, which may become less abundant 
along with their herbivorous hosts. The present study revealed how- 
ever very low parasitisms rate in both GM and non-GM maize, 
therefore they importance needs to be further evaluated. Target 
insects as the western corn rootworm and the European corn borer 
were also present. Diabrotica adults were able to persist in 
Coleoptera-resistant maize (treatments 1, 2, 7 and 8) leaves where 
its intensive predation from the ladybirds beetles were detected. Very 
low density of European corn borer was detected in all plots; its 
abundance reached a maximum of 0.2% between herbivores. 

The structure of interactions between trophic groups shows small 
if any differences between GM and non-GM maize. The complexity 
parameters together (S and L) indicate that neither Bt toxins, gly- 
phosate tolerance or the extra glyphosate treatment had any effect on 
trophic groups relative to non-GM maize hybrids. The persistence of 
a food web increases with connectance and values between 0.2 and 
0.5 suggest the existence of stable food webs'"''". Connectance in all 
GM and non-GM maize treatments was within this range of values 
with no difference between treatments (Table 3). No significant 
effects of bottom-up (F = 0.034, p = 0.96) and top-down (F = 
0.021, p = 0.98) indexes were detected between treatments. Similar 
and relatively low top-down indexes from predators to phytophag- 
ous species and from phytophagous species to plants indicate strong 
interferences. Bottom-up indexes also reveal that a strong interfer- 
ence exists between trophic groups (Table 3). 

Linear contrast tested between bottom-up and top-down effects of 
the two control maize and GM crops indicate a well fit with our initial 
prediction that similar food web exits in both non-GM and GM 
maize. A strong positive correlation between predators' top-down 
effects, herbivores top-down effects and plants bottom-upp effects in 
controls and GM maize were observed (Table 4). Similar trend were 
observed for herbivores bottom-up effects but only for treatments 7 
and 8 (Coleoptera, Lepidoptera and glyphosate resistance and 
Coleoptera, Lepidoptera and glyphosate resistance + glyphosate 
treatment) (Table 4). 

Discussion 

Several conventional studies have been previously conducted in 
commercial and experimental fields to evaluate the impact of GM 



maize crops that contain transgen from Bacillus thuringiensis (Bt 
crops) on beneficial non-target organisms'^. Field studies are useful 
tools to assess the effects of exposure on many species^'*^'', but quant- 
itative food web analyses can provide a much comprehensive 
information to understand the effect of Bt toxins to the structure 
of ecological communities. Concerns still exist that Coleopteran 
and Lepidopteran resistant GM maize may have cascading effects 
on the functioning of ecosystems and could disrupt community 
structure and ecosystem functioning^". Phytophagous and predator 
trophic groups in current study were previously considered as non- 
target organisms on which effects of Bt proteins may have possible 
direct (phytophagous species) or indirect (predators) harmful 
effects'^-'"''^'^"''''''^. These studies revealed no significant effects on 
non target organisms via direct (phytophagous) or indirect (pred- 
ator) pathway. Target pests as the western corn rootworm and the 
European corn borer were also present, nevertheless in low density. 
The extra glyphosate treatment reduced the overall weed density; 
broadleaf weeds were stiU able to persist in low mass because of 
existing seeds in soil. This reduced weed density was an important 
source for phytophagous groups that further explains the food-web 
complexity in these GM plots. Quantitative analyses of food web 
properties that includes all of the above species and more, revealed 
no negative effects of different of Bt toxins singly or in combination 
with glyphosate tolerance traits on arthropod food web structure, 
relative to non-GM maize. Evaluating effects of GM crops can be 
usefuUy examined using quantitative food web analyses as it helps to 
conduct field evaluations of effects of Bt toxins on target and non- 
target organisms as well as evaluate effects on the functional diversity 
and relations between trophic groups existing on different GM maize 
stands. However a long-term ramifications of the transgenic crops 
with multiple Bt genes and other traits are still to be examined in the 
coming decades and centuries, the present short-term analyses and 
method makes possible simple comparison between GM and non- 
GM maize stands in terms of potential shifts in vertical interactions 
of trophic groups. Both bottom-up and top-down values indicate a 
strong interference between trophic groups (Table 3). Other studies 
revealed that human modification of ecosystems can influence the 
degree of which bottom-up control is replaced by top-down con- 
troP"*. In cases when Bt crops or herbicide-resistant crops are com- 
pared with control different impact would be expected. To test this 
multiple regressions were performed between bottom-up and top- 
down values. Almost in all cases a strong positive relationship existed 
between control and GM Kbu and Ktd variables. Exceptions were 
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Figure 1 | (A-H) Food-webs of predatory and phytophagous guilds in different GM maize stands and their isogenic control. The size of each band 
represents the proportion of the trophic group. The food web for each GM and non-GM maize were made by considering the mean number of throphic 
groups/replicates/maize growing stage. Analyses were made on a subset of organisms that are the most abundant during the maize vegetation period (but 
before pollen spreading to exclude errors by cross-pollination), their cumulative abundance as trophic groups (herbivores, predators) for each GM and 
non-GM crops exceeded 95% and they were identified as test organisms in GM crops by EFSA. All food webs of GM and non-GM treatment presents the 
interactions in the same time periods. 



observed when herbivores bottom-up effects were compared be- 
tween controls and Coleopta resistant, Coleoptera and Lepidoptera 
resistant, Lepidoptera and glyphosate resistant and Lepidoptera and 



glyphosate resistanat + glyphosate treatment data (Table 4). This 
however can be explained by the lover K^u indexes (0.57) of a single 
group, Cicadellidae in these treatments (Kbu = 1.07 in the other 
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Table 3 | Food web parameter on GM and isogenic maize stands (CoSBiLab). ANOVA was used to test the effect variation on top down and 
bottom-up effects. Comparisons were made by sampling events/ treatments and their interactions for the following responses: weed 
coverage of four maize growth stages and abundance of phytophagous and predators 
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S - nr. of trophic species, L - nr. of links between tropfiic groups, C - connectance, Kbu - aritfimetic means of bottom-up indexes, Ktd - arithmetic means of top-down indexes, ns - not significant. Standard 
Errors are given in brackets 



Table 4 | Linear contrasts of bottom-up and top-down effects between the controls and GM maize crops. Bottom-up and top-down indexes 
were considered as separate variables for each trophic group 
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treatments). Lower bottom-up effects of Cicadellidae here indicates 
stronger interference and weak host dependence^^^'^^^ but these cannot 
be attributed directly to genetic modifications. Factors such as com- 
petitions between herbivores, abiotic factors (possible higher tem- 
perature, humidity, etc) may also influence herbivores bottom-up 
effects. The connectance offers complementary insight into food 
web persistence within an area, because connectance values for a 
particular trophic species reflect its foraging behaviour. As long as 
the most connected species are unlikely to go extinct, stability 
increases with connectance^^'^^^, with connectance values between 
0.2 and 0.3 are indicative of the most stable food web configurations 
across many terrestrial ecosystems^'^'^". Connectance values among 
our GM treatments and non-GM controls fell within this range 
(Table 3), suggesting the existence of stable food web on both GM 
and non-GM maize. 

Conclusions 

Analyses focused on those non-target organisms that are defined as 
test organisms in GM crops. All these organisms were in the same 
times the most abundant in both GM and non-GM maize. Target 
insects as the western corn rootworm and the European corn borer 
were also present in lower percentages. Diabrotica adults were able to 
persist in Coleopter a- resistant maize and very low density of 
European corn borer was detected, its abundance reached a max- 
imum of 0.2%. However arthropod movements between plots cannot 
be excluded, our data assessment and the realistically parameterized 
food web analyses based on a total of 243,896 arthropod individuals 
can be considered as new and useful method to assess the deleterious 
effects of Bt toxins and glyphosate on non-target arthropods. Despite 
the presence of different maize genotypes and different Bt toxins, 
glyphosate tolerance and extra glyphosate treatment, a complex and 
stable food web exist in maize ecosystems. 

Methods 

Ethics statement. For arthropod collection in maize field, no permits were required. 
Ethics approval for arthropods used for survey was not required by any specific 
committee. All animal work was conducted according to relevant national and 
international guidelines. 

Site characterization and sampling procedure. Thirty-two plots were established in 
Central Europe, near Budapest on chernozem soil in a randomized complete block 
design with each of the four GM varieties and 2 non-GM controls (Table 1) 
represented in each block. In addition, the two glyphosate tolerant varieties were 
replicated again in each block and subject to a glyphosate treatment in order to 
distinguish effects of GM genotype vs. effects of glyphosate chemical application on 
food web composition. Glyphosate in these treatments was applied in a total amount 
of 1060 g/ha in each year at the four and eight leaf stages of the maize plants. 

Each maize treatment was planted in monoculture in 625 m^ {25 m X 25 m) plots 
spaced 3 m apart within each block. Maize plants of similar maturity were planted in 
the retention zone surrounding the experimental fields to capture pollen surrounding 
the entire experimental field site to prevent pollen contamination of non-GM fields of 
the wider area around the experimental fields. Maize was planted in between late 
April and early May; and harvested between late October and early November. 

For two years, arthropod and weed were sampled weekly from April until the end of 
harvest. Because many arthropods could have moved between plots, thereby reducing 
impacts of the different treatments, samples were taken only from a 10 X 10 m area 
within the centre of each 625 m^ plot. In this way there was a distance of approxi- 
mately 18-20 meter between traps of different blocks. Three standard methods were 
used to sample arthropods. Ground -dwelling arthropods were sampled using pitfall 
traps. Three traps were placed in each plot {a total of 96 traps) and were operated 
weekly from May to October. All individuals collected were counted and assigned to a 
trophic group (herbivores and predators). Plant canopy dwelling insects were col- 
lected weekly with 30 X 20 cm Pherocon (trade mark name) yellow sticky traps. 
Again, three/plot (a total of 96), were used to collect insects from plants. Traps were 
operated for seven days in each year from May to October. In each seven-day interval, 
traps were changed and arthropods were identified in laboratory. Surveys for other 
arthropods that are not sensitive to yellow colour were assessed by visual observations 
during each sampling period. Fifteen plants per plot were randomly selected weekly 
(480 plants per assessment) and all arthropods were counted and identified to trophic 
group. Weed plant species and their area coverage were estimated in all blocks by 
visual survey on 3 X 1 m^ area of each plot as 0%, 10%, 20%, 30% ... and 100% at four 
different growth stages of maize (eight leaves stage (V8), twelve leaves stage (V12), 
vegetative stage, tasseling (VT) and reproductive stage, milk (R3)). 



Data analyses. Arthropods of similar prey preferences were assigned to the same 
trophic groups, where trophic group contains taxa that share the same set of predators 
and/or prey within a food web. This was done by direct observations of species 
interactions during the weekly plant surveys on the 480 plants (i.e. lady beetle 
predation on aphids and/or lady beetle predation on thrips etc.) and all interaction 
events noted. Designating species into trophic groups this way is a widely accepted 
convention in structural food-web studies because it reduces methodological biases 
related to uneven resolution of taxa within and among food webs^^. 

Designation of trophic links in the food web followed methods presented by^^. All 
observed interactions among the food web components were grouped according to 
the numbers of interactions (e.g., the number of times that lady beetle predation on 
aphids were observed) per sample period. To increase the sensibility of the methods 
two important factors were considered and data used for food web construction 
adjusted according to these factors: 

i. Because maize poUen moves long distances, this may resulted in the non-Bt 
plots producing some or many kernels with Bt toxins (production of Bt toxins is 
a dominant trait), the quantitative food web associated with each GM and non- 
GM maize was constructed using information on the trophic groups until 
vegetative stage, tasseling (VT). 

ii. Because presence and abundance of some or aU of the organisms (insects and 
weeds) investigated varying importantly during the vegetation period. For 
example, corn rootworms adults are absent in April, May, June, present in 
July and August, and absent again in September, October and November. 
Therefore food-webs for each GM and non-GM crops were constructed from 
data when the abundance of the most frequent species are the highest, but before 
poUen spreading (end of June, mead July). 

This resulted that the food webs were constructed by using 24,972 to 32,237 
arthropod individuals per treatment. All food webs of GM and non-GM treatment 
presents the interactions in the same time periods. 

Quantitative, weighted measures of bottom-up (Ki,u.i), top-down (Ktd.i) indexes 
and connectance (C) were computed for all trophic groups in each sampling period^^. 
Bottom-up and top-down indices emphasize vertical interactions and are used for 
analyzing food webs in which top-down and/or bottom-up processes are of particular 
interest^^. These indices emphasize vertical interactions and indicate a strong inter- 
ference between trophic groups. Large top-down effects require weak interference, 
while large bottom-up effects require both weak interference and strong prey 
dependence '^'^^i 

En 1 — \m 1 

For species "i", the first formula quantifies the bottom-up effect (Kbu.i) while the 
second formula quantifies the top-down effect (Ktdj), where n in the first formula is 
the number of predators preying species "i", d^ is the number of prey of its c* predator 
and is the bottom-up keystone index of the c^ predator. Symmetrically, m in the 
second formula is the number of prey eaten by species "i", fg is the number of 
predators of its e* prey and K^e is the top-down keystone index of the eth prey^^. 

A common Connectance (C) index was calculated which quantifies the linkage 
probability between any pair of trophic groups within food-web (L/S^) where L is the 
number of links between trophic groups and S the number of trophic groups^^. S and L 
can be derived by counting the number of trophic groups that interact and number of 
links (interactions) between them. If for example lady beetle predation on aphids were 
observed several times in a sampling period, this was considered as one interaction. 

ANOVA was used to test the effect variation on top down and bottom-up 
effects^^'^^. Bottom-up and top-down indexes were considered as separate variables 
for each trophic group. 

Linear contrasts were used to test bottom-up and top-down effects between the 
controls and GM crops. Because genetic modifications between maize plots varied 
and extra glyphosat treatments were used, impacts between bottom-up and top-down 
effects would be expected when Bt crops or herbicide -resistant crops are compared 
with control plots. Again bottom-up and top-down indexes were considered as 
separate variables for each trophic group. In this way the contrast analysis became a 
particular case of multiple regression analysis where R- square are determined as 
correlation coefficient between control and GM bottom-up and top-down values^^. 
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